close all;
c=zeros(25,1);
c(1)=1;
c(2)=3;
c(3)=0.5*3*5;
c(4)=0.5*5;
c(5)=0.5*5*7;
c(6)=0.5*3*7;
c(7)=0.125*35*9;
c(8)=0.125*30*9;
c(9)=0.125*3*9;
c(10)=0.125*63*11;
c(11)=0.125*70*11;
c(12)=0.125*15*11;
c(13)=0.0625*231*13;
c(14)=0.0625*315*13;
c(15)=0.0625*105*13;
c(16)=0.0625*5*13;
c(17)=1/16*429*15;
c(18)=1/16*693*15;
c(19)=1/16*315*15;
c(20)=1/16*35*15;
c(21)=1/128*6435*17;
c(22)=1/128*12012*17;
c(23)=1/128*6930*17;
c(24)=1/128*1260*17;
c(25)=1/128*35*17;
c=c/4/pi;

d=zeros(21,1);
d(1)=0;
d(2)=3;
d(3)=3*5;
d(4)=0.5*15*7;
d(5)=0.5*3*7;
d(6)=0.125*140*9;
d(7)=0.125*60*9;
d(8)=0.125*315*11;
d(9)=0.125*210*11;
d(10)=0.125*15*11;
d(11)=0.0625*1386*13;
d(12)=0.0625*1260*13;
d(13)=0.0625*210*13;
d(14)=1/16*429*7*15;
d(15)=1/16*693*5*15;
d(16)=1/16*315*3*15;
d(17)=1/16*35*1*15;
d(18)=1/128*6435*8*17;
d(19)=1/128*12012*6*17;
d(20)=1/128*6930*4*17;
d(21)=1/128*1260*2*17;
d=d/4/pi;

format long;
disp(c);
disp(d);

figure;
semilogy(c,'o-','linewidth',2);
figure;
semilogy(d,'o-','linewidth',2);